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Abstract 

Background: In addition to alterations concerning the expression of oncogenes and onco-suppressors, melanoma 
is characterized by the presence of distinctive gangliosides (sialic acid carrying glycosphingolipids). Gangliosides 
strongly control cell surface dynamics and signaling; therefore, it could be assumed that these alterations are linked 
to modifications of cell behavior acquired by the tumor. On these bases, this work investigated the correlations 
between melanoma cell ganglioside metabolism profiles and the biological features of the tumor and the survival 
of patients. 

Methods: Melanoma cell lines were established from surgical specimens of AJCC stage III and IV melanoma 
patients. Sphingolipid analysis was carried out on melanoma cell lines and melanocytes through cell metabolic 
labeling employing [3-^H]sphingosine and by FACS. N-glycolyl CMS was identified employing the 14 F7 antibody. 
Gene expression was assayed by Real Time PGR. Gell invasiveness was assayed through a Matrigel invasion assay; 
cell proliferation was determined through the soft agar assay, MTT, and pH] thymidine incorporation. Statistical 
analysis was performed using XLSTAT software for melanoma hierarchical clustering based on ganglioside profile, 
the Kaplan-Meier method, the log-rank (Mantel-Gox) test, and the Mantel-Haenszel test for survival analysis. 

Results: Based on the ganglioside profiles, through a hierarchical clustering, we classified melanoma cells isolated 
from patients into three clusters: 1) cluster 1, characterized by high content of GM3, mainly in the form of N-glycolyl 
GM3, and GD3; 2) cluster 2, characterized by the appearance of complex gangliosides and by a low content of GM3; 
3) cluster 3, which showed an intermediate phenotype between cluster 1 and cluster 3. Moreover, our data demonstrated 
that: a) a correlation could be traced between patients' survival and clusters based on ganglioside profiles, with cluster 
1 showing the worst survival; b) the expression of several enzymes (sialidase NEU3, GM2 and GMl synthases) involved 
in ganglioside metabolism was associated with patients' survival; c) melanoma clusters showed different malignant 
features such as growth in soft agar, invasiveness, expression of anti-apoptotic proteins. 

Conclusions: Ganglioside profile and metabolism is strictly interconnected with melanoma aggressiveness. Therefore, 
the profiling of melanoma gangliosides and enzymes involved in their metabolism could represent a useful prognostic 
and diagnostic tool. 
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Background 

Melanoma is the most lethal form of skin cancer and ac- 
counts for at least 48,000 deaths worldwide annually 
[1,2]. Clinical outcome of melanoma largely depends on 
the tumor stage upon first diagnosis. Currently, melan- 
oma staging is based on the guidelines published by the 
American Joint Committee on Cancer (AJCC) in 2009 
[3] and advises the employment of histopathological and 
clinical criteria. Nevertheless, this system is limited in its 
ability to provide a precise prognosis: a large number of 
patients with similar or identical clinical and histo- 
pathological features has different clinical outcome, from 
being cured to death [4]. Also, the identification of some 
mutations concerning oncogenes, including BRAF, 
NRAS, KIT, has been shown to be very useful to predict 
the response to therapy, but it is not clear whether the 
determination of the mutational status could offer any 
prognostic value [5]. Many attempts have been done to 
identify molecular markers or gene and protein signa- 
tures predicting clinical outcome in melanoma but, so 
far, none of them has been sufficiently validated [6-9] . 

Among the distinctive molecular markers displayed by 
melanoma cells, many studies identified specific ganglio- 
sides not detectable in normal melanocytes [10]. A pivotal 
role played by gangliosides was recognized in the develop- 
ment of melanoma and in its biological features [11]. In 
particular, GD3 and derivatives [12,13], which are also 
present at high levels in proliferating cells [14-16], have 
been demonstrated to promote melanoma cell proliferation 
and invasion [17]. Moreover, de-N-acetyl GM3 (d-GM3) (a 
variant of GM3 with a free amino group at the 5 position 
of sialic acid instead of the acetyl group) was found in 
melanoma and it was demonstrated that it enhanced 
cell migration and invasion [18]. Also, N glycolyl GM3 
(Neu5Gc-GM3) (a variant of GM3 that contains N- 
glycolylneuraminic acid instead of N-acetylneuraminic acid) 
was recognized in melanoma, even its role in this disease is 
still obscure [19]. The alteration of the expression/catalytic 
activity of glycosyltransferases and enzymes involved in 
gangliosides metabolism could be critical factors in deter- 
mining the melanoma cell ganglioside composition. In par- 
ticular, the activation of GD3 synthase is linked to the 
increment of GD3 in melanoma cells [11,20]. Similarly, the 
plasma membrane sialidase NEU3 has been shown to be 
highly expressed in human melanoma cell lines [21]. 

Prompted by these data, we sought to investigate the 
ganglioside metabolism profile of metastatic melanoma 
cell lines established from patients. Our results demon- 
strated that: a) melanomas displayed different ganglio- 
side patterns and three clusters of tumors could be 
identified; b) a correlation could be traced between pa- 
tients' survival and melanoma ganglioside profiles; c) the 
expression of several enzymes involved in ganglioside 
metabolism was associated with patients' survival; 



d) melanoma clusters identified on the basis of ganglio- 
side profile exhibited different features determining mel- 
anoma malignancy. 

Methods 

Cell cultures 

Melanoma cell lines were established from surgical spec- 
imens of AJCC stage III and IV melanoma patients 
admitted to Fondazione IRCCS Istituto Nazionale 
dei Tumori, Milan [22,23]. Molecular and biological 
characterization of the cell lines has been reported pre- 
viously [24]. All cell lines were maintained as described 
[25]. All patients were informed about the scope and 
methods and delivered a written informed consent for 
the use of the surgical samples to establish cell lines. 
The study was approved by the Ethics Committee of the 
University of Milan and was performed according to the 
Declaration of Helsinki. 

Clones 2/14 and 2/21 were isolated from a single human 
metastatic melanoma cell line, as described [26,27]. 
NHEM-Ad and NHEM-Neo were purchased by Lonza 
(Basel, Switzerland) and PromoCell (Heidelberg, Germany), 
and maintained in mMGM-4 medium (Lonza). 

Sphlngolipid analysis 

Sphingolipid analysis was carried out through cell me- 
tabolic labeling with [3-^H]sphingosine (PerkinElmer, 
Waltham, MA, USA) [28]. In order to assay the hypoth- 
esis that NeuSGc-glicans could be incorporated from the 
culture medium and then employed for the synthesis of 
GM3, before [3-^H]sphin gosine labeling, melanoma L6 
cells were pre-incubated in the reduced-serum medium 
OptiMEM (Life Technology, Carlsbad, CA, USA) for 
5 days. 

Ganglioside and neutral sphingolipid extracts were an- 
alyzed by HPTLC carried out with the solvent systems 
chloroform/methanol/0.2% CaCl2 55:45:6 (v/v) and 
chloroform/methanol/water 110:40:6 (v/v), respectively. 
To separate Neu5Gc-GM3 from Neu5Ac-GM3, HPTLC 
was carried out using the solvent system chloroform/ 
methanol/0.2% CaCl2/5 N NH3 50:42:6:4 (v/v). The 
sphingolipid pattern was determined and quantified by 
radiochromatoscanning (Betaimager 2000, Biospace, Paris, 
France) [28,29]. 

Endogenous sphingolipid analysis performed to stan- 
dardize metabolic labeling was performed as previously 
described [30]. 

Ganglioside standards were kindly given by Prof. 
Sonnino, University of Milan. 

Immunostaining of HPTLC 

After the chromatographic separation of gangliosides, 
the plates were soaked in acetone plus 0.1% polyisobu- 
tylmethacrylate. After drying and blocking with PBS -4% 
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milk, the plates were incubated with 5 (ig/ml of anti- 
Neu5Gc-GM3 murine 14 F7 antibody [31,32], overnight, 
and, then, with a horseradish peroxidase-conjugated 
anti-mouse IgG antibody (Santa Cruz Biotechnology, 
Santa Cruz, CA, USA). The reaction was stained with 
the SuperSignal West Pico Chemiluminescent Substrate 
(Thermo Fisher Scientific, Dallas, MA, USA). 

Immunohistochemical analysis 

Cytospin preparations of melanoma cell lines were fixed 
with formalin for 10 min at room temperature and then 
incubated with 14 F7 mAb, as described [19]. The slides 
were then incubated with goat anti-mouse biotinylated 
secondary antibody, followed by streptavidin/HPR (Dako, 
Glostrup, Denmark), for 30 min at each step. The enzyme 
activity was detected with a commercial solution of 3, 3'- 
Diaminobenzidine ( Sigma- Aldrich, St Louis, MO, USA). 
Cells were counterstained with Mayers Hematoxylin 
(Bio-Optica, Milan, Italy). Slides subjected to the same 
treatments, but without incubation with primary mAb, 
were used as negative controls. Images were acquired 
at X 20 with an Axiovert 100 microscope (Zeiss, Carl 
Zeiss, Thornwood, NY) equipped with a digital camera 
(AxioCam MrC5, Zeiss). 

Phenotypic profile of melanoma cell lines 

Cell surface expression of GM3, GD2 and GD3 was 
determined in representative melanoma cell lines be- 
longing to different clusters, by staining with the follow- 
ing mAbs: anti GM3 (Neu5Ac-GM3) M2590 (Cosmo 
Bio Co, LTD, Tokyo, Japan), anti-human disialoganglio- 
side GD2 (BD Pharmingen, San Diego, CA, USA), anti- 
ganglioside GD3 (R24) (Abeam Inc., Cambridge, UK), 
followed by a FITC-conjugated F(ab')2 fragment goat anti- 
mouse Ab (Jackson ImmunoResearch, West Grove, PA) as 
described [33]. Samples were acquired by a fluorescence- 
activated cell sorting (FACS)-Calibur cytofluorimeter (BD 
Biosciences) and analyzed by the Flowjo software (Tree 
Star, Ashland, OR, USA). Results were expressed as per- 
centage of positive cells and as mean fluorescence intensity 
after subtraction of mean fluorescence intensity in cells 
stained with secondary antibody only. 

In order to compare the phenotype of different tumors, 
all main FACScan operative settings were kept constant 
throughout the whole study. 

PGR and Real-time RT-PCR 

CMAH expression was assayed in cDNA isolated from 
melanocytes and melanoma cells through RT-PCR. The 
genomic region corresponding to the 92-bp deletion in 
the human CMAH gene was amplified by PCR using 
genomic DNA extracted from melanoma cells and from 
murine C2C12 ceUs, as positive control. Real Time PCR 



was performed as previously reported [34]. Primer se- 
quences are reported in Table 1. 

Invasiveness assay 

Melanoma cell invasiveness was assayed through a 
Matrigel assay, as previously reported [35]. 

Proliferation assays 

Soft agar, thymidine incorporation, and MTT assays 
were performed as previously reported [34]. 

Sialidase activity toward endogenous cell gangliosides 

To test the abflity of melanoma plasma membrane si- 
alidase to act on endogenous cell gangliosides, melan- 
oma cell particulate fractions, obtained as previously 
described [36], were incubated with [3-^H]sphingosine 
labeled ceU gangliosides at pH 3.8, overnight, at 37°C. 
Then, gangliosides were fractionated by HPTLC and vi- 
sualized by radiochromatoscanning. 

Statistical analysis 

Hierarchical clustering was done using XLSTAT soft- 
ware. Survival analysis was carried out by the Kaplan- 
Meier method; survival curves were compared by the 
log-rank (Mantel- Cox) test; hazard ratios were com- 
puted by the Mantel-Haenszel test. Ganglioside compos- 
ition. Real Time PCR data, [^H] thymidine incorporation 
were compared using Student t-test and 1-way ANOVA. 

Results 

Melanoma cells can be clustered on the basis of their 
ganglioside profile 

Sphingolipid pattern was determined in 23 short-term 
ceU lines established from 5 primary VPG (vertical 
growth phase) melanomas and from 18 lymph nodes 
metastases as well as in adult (NHEM-Ad) and neonatal 
(NHEM-Neo) melanocytes. To this end, [3-^H]sphingo- 
sine was administrated to the cells leading to an exten- 
sive labeling of all sphingolipids [37]. The treatment 
was performed in order to reach a metabolic steady 
state that corresponded to endogenous non-radiolabelled 
cell sphingolipid pattern. NHEM-Ad and NHEM-Neo 
melanocytes ganglioside profile was mainly composed of 
GM3 and GD3 (in traces in NHEM-Ad cells and more 
abundant in NHEM-Neo cells) with traces of sialyl- 
paragloboside (SPG) (Additional file 1). 

The distinct ganglioside profiles of melanoma cells 
allowed us to classify them into three main groups 
(Figure lA). The ganglioside HPTLC pattern of NHEM- 
Ad and NHEM-Neo melanocytes and of representative 
cell lines of each cluster is shown in Figure IB. Melanoma 
cells belonging to cluster 1 showed a strong similarity with 
neonatal melanocytes, with GM3 as the main ganglioside 
and an increased amount of GD3 (Additional file 1, 
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Table 1 Primers used for gene expression 



Gene 


Forward primer 


Reverse primer 


GD3 synthase {STSSIA!) 


5'-CAGCATAATOGGCAAAGGT-3' 


5'-ATOGCACCAACATCTGACA-3' 


NEUl 


5'-CCTGGATATOGCACTGAA-3' 


5'-CATCGCTGAGGAGACAGAAG-3' 


NEU2 


5'-AGAAGGATGAGCACGCAGA-3' 


5'-GGATGGCAATGAAGAAGAGG-3' 


NEU3 


5'-GGCTOmGGGTGmG^-3' 


5'-CATCGCTGAGGAGACAGAAG-3' 


NEU4 


5'-ACCGCCGAGAGIGI 1 1 IGG-3' 


5'-CGTGGTCATCGCTGTAGAAGG-3' 


GM3 synthase {ST3GAL5) 


5'-CCCTGAACCAGTOGATG^-3' 


5'-CATOCTOAAGCCAGTOA-3' 


GMl synthase {B3GAW4) 


5'-CGCCTOCAGGACTCCTACC-3' 


5'-CCGTCTOAGGACGTATCGG-3' 


GM2 synthase {B4GALNT1) 


5'-TCTCACTGGAGAGGGTCAGG-3' 


5'-CGGACTGTGTCTGCTGTG^-3' 


GDI a synthase {ST6GALNAC4) 


5'-TCTACCACCCAGCCTOATC-3' 


5'-TAGTGGTGCCAGTOCCm-3' 


STSGolVI 


5'-TOCCCTATGGGATGAGAAC-3' 


5'-CCTCCA™CCAACCACCAC-3' 


ST3GalNac6 


5'-CGCCGGAGAGAAATGAGTA-3' 


5'-CCACTOTOAGGTOACAGG-3' 


SPl 


5'-GGCTACCCCTACCTCAAAGG-3' 


5'-CACAACATACTGCCCACCAG-3' 


CMAH 


5 '-G A AmATOTCCCmGC-3 ' 


5'-TATGATGCCCTOCTGTC-3' 


CMAH deletion 


5'-AGCmCCGGAGGGTOATA-3' 


5'-TGTCAGATGCACAAGCACAA-3' 


ACTB (P actin) 


5'-CGACAGGATGCAGAAGGAG-3' 


5'-ACATCTGCTGGAAGGTGGA-3' 



Figure IB), Instead, ganglioside profiles of cells belonging 
to clusters 2 and 3 were significantly different from that of 
melanocytes, revealing the appearance of novel species 
such as GM2, GDla, GD2, GTlb in cluster 2, and GM2, 
GDla, GD2 in cluster 3, in addition to GM3 and GD3 
(Additional file 1, Figure IB). No significant differences 
were found considering primary lesions versus metastases 



or the clinical stage of the disease at the time when 
tumors were surgically removed and cell lines established 
(AJCC stage III versus stage IV) (Additional file 1). 

Despite these huge differences regarding the ganglioside 
profiles of the three clusters, neutral sphingolipid pattern 
was quite similar also to adult and neonatal melanocytes 
pattern, except for a 30-40% decrease of ceramide amount 
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Figure 1 Association of melanoma cell ganglioside profile with patients' survival. (A) Hierarchical clustering analysis of melanoma cell lines on 
the basis of ganglioside profiles. Dotted line represents the level of significant dissimilarity. (B) HPTLC separation of gangliosides and (C) of neutral 
sphingolipids of three representative melanoma cell lines (L6, L34, L4) belonging to clusters 1,2, and 3, in comparison to adult and neonatal 
melanocytes. (D) Kaplan-Meier survival analysis of patients according to the three clusters identified. Survival was evaluated as time from first diagnosis. 
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in cells classified in cluster 1 in comparison to cells classi- 
fied in cluster 2 (P < 0.05) (Figure IC). 

Kaplan-Meier survival analysis and log-rank test demon- 
strated that patients whose tumors fell in cluster 1 had a 
significantly worse overall survival, measured from the 
time of the first diagnosis, in comparison to patients 
whose tumors were classified in cluster 3 and, above all, in 
cluster 2 (P < 0.05) (Figure ID). Patients belonging to clus- 
ter 1 had a hazard ratio of 4.23 in comparison to patients 
belonging to cluster 2 (P < 0.05), and a hazard ratio of 2.91 
as compared with patients belonging to cluster 3 (P < 0.05). 

The ganglioside phenotypic profile of representative cell 
lines for each cluster is reported in Figure 2. A low hetero- 
geneity was observed concerning the content of GD3 and 
GD2 in some cell lines. 

A high amount of N-glycolyl GM3 characterizes melanoma 
cluster 1 cells 

We recorded a singular trend showing a marked decrease 
of the average GM3 per cent content from melanoma 
cluster 1 cells to cluster 3 (-58.3%; cluster 3 versus cluster 
1, P < 0.001), and to clusters 2 (-76.5%; cluster 2 versus 
cluster 1, P < 0.001) (Figure 3A). This trend was associated 
with the decrease of GM3 synthase expression in clusters 
2 and 3, in comparison to cluster 1 (-85.7%, cluster 2 ver- 
sus cluster 1; -82.1% cluster 3 versus cluster 1; P < 0.001) 
(Figure 3B). To confirm the key role played by GM3 in 
metastatic melanoma behavior, we analyzed the ganglio- 
side composition of two melanoma clones (2/14 and 2/21) 
isolated from a single melanoma cell line [26]. These two 
clones had markedly different molecular and biological 
features [26,27], including adhesion molecule phenotype, 
degree of differentiation, recognition by antigen-specific T 
cells, presence of mutated NRAS (clone 2/14) or mutated 
BRAF (clone 2/21), growth in immunodeficient mice (en- 
hanced in clone 2/14) [26,27]. Significantly, the ganglioside 
profile of clone 2/14 consisted only of GM3 (Figure 4A). 

We assessed that a high amount of GM3 was present in 
clone 2/14, as well as in melanoma cell lines, as Neu5Gc- 
GM3. Neu5Gc-GM3 was detected also in melanocytes. 
d-GM3 was not present (Figure 4A). The increase of GM3 
in melanoma cluster 1 cells reflected an increased content 
of Neu5Gc-GM3 (90% of total cellular GM3) in compari- 
son to other clusters (Figure 4B). The identification of 
Neu5Gc-GM3 was validated by the employment of the 
Neu5Gc-GM3-specific antibody 14 F7 [19] in HPTLC 
immunostaining (Figure 4C) and in immunohistochemical 
analysis, which confirmed its enhanced amount in cluster 1 
melanoma cells (Figure 4D). Because it was reported that 
N-glycolylneuraminic acid synthesis could not occur in 
human cells due to a 92 bp sequence deletion in the 
CMAH gene that blocks the synthesis of a fijinctional form 
of CMP-NeuAc hydroxylase (CMAH), the enzyme respon- 
sible for the conversion of CMP-N-acetylneuraminic acid 



into CMP-N-glycolylneuraminic acid [38,39], firstly we ex- 
plored if this genetic alteration was reverted in melanoma 
cells. We recorded CMAH expression as mRNA in all 
melanoma cell lines, in contrast to melanocytes where it 
was null (Figure 5A). Nevertheless, we demonstrated that 
the CMAH gene maintained the 92 bp sequence deletion 
typical of human ceUs [40], also in melanoma ceUs 
(Figure 5B); therefore, on the basis of previous evidences 
[41], this gene cannot be translated into an enzymatically 
active protein. Many reports demonstrated that some can- 
cer cells become able to incorporate dietary CMP-N- 
acetylneuraminic acid at a faster rate from biological fluids 
and therefore, they incorporate it into gangliosides. Effect- 
ively, we demonstrated that cultivating melanoma cells in a 
serum-reduced medium for 5 days, the content of 
Neu5Gc-GM3 was reduced by 25% (P < 0.01) (Figure 5C). 
Therefore, we concluded that melanoma cells do not gain 
the capability to synthetize CMP-N-acetylneuraminic acid 
but they can acquire it from external molecules. 

Correlation between glycosyltransferases expression and 
patients' survival 

We confirmed an overall higher expression of GD3 
synthase (11.5-1.5-fold) and plasma membrane sialidase 
NEU3 (7-1.5-fold) in all melanoma cells than melanocytes 
(Figure 6A,B) [11,21]. The increase of the expression of 
these two enzymes was reasonably induced by the 
enhanced expression of the transcriptional factors Spl 
(average increase: 100%, melanoma cells versus adult mela- 
nocytes; P < 0.001) (Figure 6C), as it was reported that both 
NEU3 [42] and GD3 synthase gene promoters [43] are 
regulated by them. NEU3 was the unique sialidase to be 
over-expressed by melanoma cells; lysosomal sialidase 
NEUl expression varied inconsistently among melanoma 
cell lines and melanocytes; cytosolic sialidase NEU2 was 
not expressed in both melanocytes and melanoma cells; 
sialidase NEU4 appeared to be down-regulated in all 
melanoma cells except LI, L34, L39, L40, in comparison to 
melanocytes. As melanoma cells expressed high levels of 
GM3, the increase of NEU3 expression appeared a paradox. 
Therefore, we tried to assess if endogenous GM3 could be 
accessible to NEU3 activity in melanoma cells. To this end, 
after metabolic labelling with [3-^H]sphingosine, gang- 
liosides were extracted from a representative melanoma cell 
line belonging to cluster 1 (L6) and incubated with mem- 
branes obtained by the same cell line, in order to perform 
an exhaustive treatment. We demonstrated that NEU3 was 
able to hydrolyze and convert into lactosylceramide only 
41% of plasma membrane GM3: the aliquot of Neu5Gc- 
GM3 was particularly resistant to sialidase action (-25- 
30%), while N acetyl GM3 (Neu5Ac-GM3) was entirely 
hydrolyzed. Also, GD3 and SPG decreased by 22% 
and 34%, respectively, after incubation (pre-incubation 
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Figure 2 Ganglioside phenotypic profile of melanoma cell lines. Cell surface staining of GM3, GD2, and GD3 by mAbs: anti Neu5Ac-GM3 (M2590); 
anti-GD2, anti-GD3 (R24), determined in representative melanoma cell lines belonging to different clusters (L40 and L23 for cluster 1; LI and L34 for 
cluster 2; L2 and LI 6 for cluster 3). 



gangliosides versus post-incubation gangliosides; P < 0.01) 
(Figure 6D,E). 

We recorded in some melanoma cell lines grouped, par- 
ticularly, in clusters 2 and 3, the expression of GMl and 



GM2 synthases that were not expressed by melanocytes 
(Figure 7A). The expression of GMl and GM2 synthases 
was related to the appearance of GMl and GM2 among 
the gangliosides of these cells (Additional file 1). Kaplan- 
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Figure 3 Correlation between GIVIS amount and melanoma cell malignancy. (A) GM3 content and (B) Real Time PGR analysis of the mRNA 
expression of GM3 synthase in melanoma cells and in melanocytes. Values are the mean + S.D. of four independent experiments. 



Meier survival analysis and log-rank test demonstrated that 
patients whose tumors showed a relative mRNA expression 
of GMl synthase, normalized on |3 actin, higher than 1 had 
a significantly better overall survival, measured from the 
time of the first diagnosis, in comparison to patients whose 



tumors showed a less relative mRNA expression of this 
enzyme (P < 0.01) (Figure 7B). Patients carrying tumors 
endowed with GMl synthase low expression had a hazard 
ratio of 4.499 in comparison to other patients (P < 0.01). 
Similarly, patients whose tumors showed a relative mRNA 
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Figure 4 Correlation between Neu5Gc-GM3 and melanoma cell malignancy. (A) HPTLG separation of gangliosides of two cell clones (2/14 
and 2/21) isolated from a single patient. (B) HPTLG separation of Neu5Ac-GM3 and Neu5Gc-GM3 in NHEM-Ad, NHEM-Neo, clone 2/14, and in 
melanoma cell lines representative of cluster 1 (L6), cluster 2 (L34), cluster 3 (L4). (C) 14 F7 antibody immunostaining of gangliosides separated by 
HPTLG (NHEM-Ad, NHEM-Neo, L6 representative of cluster 1, L34 representative of cluster 2, L4 representative of cluster 3. (D) Immunohistochemical 
analysis of cytospin preparations of L6 representative of cluster 1 and L34 representative of cluster 2, stained with 14 F7 antibody. Original 
magnification, x20 (Axiovert 100). Similar results were obtained for other melanoma cell lines according to the clustering. 
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expression of GM2 synthase, normalized on |3 actin, 
higher than 1.5 had a significantly better overall survival, 
measured from the time of the first diagnosis, in compari- 
son to other patients (P < 0.05) (Figure 7C). Low or null 
expression of GM2 synthase entailed for patients a hazard 
ratio of 3.5 in comparison to high expression (P < 0.01). 

Overall, melanomas that showed high expression of both 
GMl and GM2 synthases were associated with an overall 
survival longer than 5 years, while the lack of both these 
glycosyltransferases increased the hazard ratio to 6.325 
(P < 0.01) (Figure 7D). All melanoma cluster 2 cells and the 
majority of melanoma cluster 3 cells expressed both GMl 
and GM2 synthases (Figure 7A). 

We also analyzed the expression of other glycosyl- 
transferases involved in ganglioside metabolism such 
as GDla synthase, ST3GalVI, ST6GalNAc6 (data not 
shown) but we did not identify any correlation with 
patients' survival. 



Correlation between melanoma cell ganglioside profile 
and growth/survival and invasive potential 

The anchorage-independent growth potential of melan- 
oma cells was assessed through a soft agar assay. Melan- 
oma 1 cells were able to form more and larger colonies in 
soft agar than were melanoma cluster 2 cells (4.1-fold; P < 
0.01), and than were melanoma cluster 3 cells (2-fold; P < 
0.01) (average colony forming efficiency: cluster 1: 53%, 
cluster 2: 13%; cluster 3: 27%) (Figure 8A). This different 
potential of growth in soft agar did not reflect a similar in- 
hibition of proliferation capability. In fact, [^H] thymidine 
incorporation was similar between melanoma cluster 1 
and 2 cells and even higher in cluster 3 (Figure 8B); MTT 
assay did not reveal any difference among the three clus- 
ters (Figure 8C). Therefore, we could not conclude that 
the ganglioside profile of the three clusters was determin- 
ant for the proliferation rate but only for the tumorigen- 
icity measured as capability to grow in soft agar. 

In the Matrigel invasion assay, cluster 1 cells revealed to 
be more invasive than cluster 2 and 3 cells (+37%; P < 0.01) 
(Figure 8D). 

The resistance to apoptosis evidenced by the expression 
of BCL-2 and BCL-XL was clearly different between the 
three clusters: melanoma cells grouped in cluster 1 
showed a higher expression of BCL-2 than melanocytes 
and also than melanoma cells grouped in clusters 2 and 3 
(+150% versus adult melanocytes, P < 0.001; +150% versus 
cluster 2, P < 0.001; +900% versus cluster 3, P < 0.001) and 
BCL-XL (+38% versus adult melanocytes, P < 0.01; +57% 
versus cluster 2, P < 0.001; +38% versus cluster 3, P < 0.01) 
(Figure 8E). 

Discussion 

The expression of atypical gangliosides deeply marks 
melanoma [20]. This study demonstrated that melano- 
mas could be classified based on their ganglioside profile 
and that there is a significant correlation between the 
ganglioside pattern displayed by melanomas and their 
biological and clinical behavior. Ganglioside biosynthesis 
occurs via two major pathways designed "a" (GM2, 
GMl, GDla, and GTla) and "b" (GD3, GD2, GDlb, and 
GTlb) with a common precursor (GM3) (Figure 9). In- 
triguingly, melanoma ceUs that we grouped in cluster 1, 
retrieved and exacerbated ganglioside biosynthetic path- 
ways typical of neonatal melanocytes, characterized by a 
high content of GD3 ("b" series) and GM3. GD3 could be 
considered a typical fetal ganglioside [44]. In melanoma 
cells, GD3 was demonstrated to enhance the adhesion and 
proliferation properties by recruiting integrins and signal 
molecules to lipid rafts [20,45]. Furthermore, the 9-0 
acetylation of GD3 (to form 9-OAc-GD3) was described 
in melanoma [46]: this chemical modification suppresses 
the pro-apoptotic flanction of GD3 [47]. Moreover, 9-OAc- 
GD3 is an important regulatory molecule being involved in 
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signal transduction, regulation of growth, and differenti- 
ation [47]. Regarding GM3, it should be noted that, in 
cluster 1, it was principally in the form of Neu5Gc-GM3 
instead of Neu5Ac-GM3. N-glycolylneuraminic acid syn- 
thesis has been considered to be absent in normal human 
tissues due to a partial deletion of the gene coding for 
CMAH [38,39]. We demonstrated that melanoma cells 
characteristically expressed CMAH cDNA. Nevertheless, 
also in melanoma cells, CMAH gene carried the 92 bp 
deletion [40]. Thus, it cannot be hypothesized that the 



synthesis of N-glycolylneuraminic acid could be rescued 
in melanoma cells. Otherwise, we determined that 
Neu5Gc-GM3 amount in melanoma cells depended on 
the external medium. Therefore, it can be supposed that 
melanoma cells, in particular those grouped into cluster 1, 
could incorporated N-glycolylneuraminic acid in a faster 
way, as described for other tumors [38]. It has been 
previously showed a correlation between CMAH mRNA 
expression and an increased ability of human stem cells to 
bind and possibly ingest exogenous N-glycolylneuraminic 
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acid [48]. Therefore, it could be speculated that the 
appearance of CMAH protein, even though enzymati- 
cally inactive, could play a role in the incorporation of 
external N-glycolylneuraminic acid. The incorporation of 
N-glycolylneuraminic acid in vivo could possibly occur 
from biological fluids. In fact, the presence of NeuSGc- 
GM3 in melanoma surgical specimens was previously 
largely documented [19,32,49]. Moreover, melanoma 
cluster 1 cells were characterized by the highest levels 
of expression of GM3 synthase: therefore, also N- 
glycolylneuraminic acid could be recruited faster into 
GM3, in these conditions. The increase of NeuSGc- 
GM3 was demonstrated to enhance the aggressive 
behavior of cancer. Mouse B16 melanoma cells loaded 
with Neu5Gc-GM3 showed more lung metastases and 
increased tumor burden in subcutaneous implants [50]. 
Moreover, on the contrary of Neu5Ac-GM3, NeuSGc- 
GM3 did not inhibit EGF-induced EGFR phosphorylation 
in A431 cells [51] and down-modulated CD4 molecule 
in T lymphocytes [52] and dendritic cell differentiation 
[53]. To further confirm the involvement of GM3 and 
Neu5Gc-GM3 in melanoma aggressiveness, we demon- 
strated that a clone isolated from a single melanoma cell 



line and endowed with higher adhesive and invasive 
potential, showed a gangUoside pattern composed only 
by GM3, mainly in the form of Neu5Gc-GM3. Neu5Gc- 
GM3 was not efficiently hydrolyzed by plasma mem- 
brane sialidase and this explains why melanoma cells 
could accumulate it despite the high expression of 
NEU3, In our melanoma cells, we could not identify d- 
GM3 in contrast to Liu et al. who quantified d-GM3 as the 
major GM3 species present in melanoma [18]. Summing 
up, melanoma cluster 1 cells showed a ganglioside profile 
that appeared to confer the highest malignant properties in 
terms of growth in soft agar medium, in vitro invasiveness, 
and expression of anti-apoptotic proteins. Also, the increase 
of BCL-2 and BCL-XL could be the direct consequence of 
an altered control of growth factor receptors such as EGFR, 
due to the increase of both GD3 [54] and Neu5Gc-GM3 
[51] on the plasma membrane. 

Melanoma cluster 2 ceUs completely revolutionized 
the usual pathways of ganglioside biosynthesis present in 
melanocytes and acquired the expression of GMl and 
GM2 synthases. Their ganglioside pattern was enriched 
by the presence of complex gangliosides such as GDla 
and GTlb, belonging to the "a" pathway, which were not 
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Figure 8 Growth, invasiveness, and survival potential of melanoma clusters. (A) Images of colonies formed by representative melanoma 
cells grouped in cluster 1 (L6), 2 (LI), and 3 (L4) in the soft-agar assay after 3 weeks of incubation at 37°C Original magnification x4 (Olympus 
1x50). Three replicate experiments were performed. Similar results were obtained for other melanoma cell lines according to the clustering. (B) 
Average value of proliferation of melanoma cells grouped in cluster 1, 2, and 3 assessed using pH] thymidine incorporation and (C) assay. 
NS: non significant. Four replicate experiments were performed. (D) Average invasiveness potential assessed by Matrigel invasion assay of melanoma 
cells grouped in cluster 1, 2, and 3. Three replicate experiments were performed. (E) Average mRNA expression of BCL-2 and BCL-XL in melanoma cell 
lines classified in clusters 1, 2, and 3 and in adult melanocytes. Dotted line represents BCL-2 and BCL-XL expression level in adult melanocytes. Values 
are the mean ±S.D. of four independent experiments. 
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revealed in melanocytes, and to the "b" pathway. It was 
previously demonstrated that the transfection of GMl 
synthase cDNA into the melanoma cell line SK-MEL-37 
gave rise to the neo-expression of GDlb, GTlb, and 
GMl and, in parallel, reduced cell growth and invasion 
[55]. In melanoma cluster 2 cells, the expression of GM3 
synthase was similar to adult melanocytes and, also, the 
content of GM3 and Neu5Gc-GM3 was significantly less 
than melanoma cluster 1 cells. Melanomas classified in 
cluster 3 represented a middle situation between that 
displayed by melanoma cluster 1 and melanoma cluster 
2 with low Neu5Gc-GM3 levels and GM3 synthase ex- 
pression and the presence of GDla. Accordingly to these 
different ganglioside profiles, also the in vitro malignant 
behavior showed by these two clusters was clearly weak- 
ened. Thus, it appears realistic that the synthesis of dif- 
ferent gangliosides, altering the recruiting and activation 
of signal molecules, could reflect on the biological prop- 
erties of melanoma. Significantly, also the clinical behav- 
ior was very different: we traced a correlation between 
melanoma clusters and patients' survival. In particular, 
the synthesis of GM3/GD3 and of "a" complex gangliosides 
appeared to have a worst and good prognostic value, re- 
spectively. Intriguingly, melanomas that retrieved a ganglio- 
side metabolism typical of immature melanocytes (cluster 
1), showed the worst outcome (median survival lower than 
10 months). It could be supposed that this phenotype 
corresponds to a higher motile behavior: in fact, neural 
crest-derived melanoblasts are highly migratory cells [56]. 
Instead, complex gangliosides, such as GDla, usually de- 
note more differentiated cells [57]. In the past, pioneering 
studies have indicated that the GM3: GD3 ratio could be 
followed for the prognosis and therapeutic management of 
melanoma [58]. Our results demonstrated that the overall 
cell ganglioside profile and metabolism could be more 
important than the presence of a single ganglioside type lil<e 
GD3 or N-glycolyl GM3 for cell biologic features and, 
therefore, for prognosis. 

Sialidase NEU3 and GD3 synthase genes were signifi- 
cantly up-regulated in melanomas in comparison to 
melanocytes, possibly as a direct consequence of the 
increase expression of the transcriptional factor Spl, but 
we cannot record any significant correlation with 
patients' survival. Sialidase NEU3 could be involved in 
melanoma malignancy decreasing the levels of Neu5Ac- 
GM3. It should also be noted that previous papers dem- 
onstrated that NEU3 is able to modulate cell signaling 
also inducing minimal changes in ganglioside profile 
[21,59] and that NEU3 up-regulation could be counter- 
balanced by other modifications concerning the expres- 
sion of enzymes involved in ganglioside metabolism, 
including GM3 synthase, inducing not expected effects 
on gangliosides [60]. Therefore, the impact of NEU3 up- 
regulation in melanoma could be complex and should 



be further investigated in forthcoming works. In sum- 
mary, the following factors appeared to be determinant 
for melanoma behavior and could have a significant 
prognostic value: a) the levels of GM3 and GD3; b) the 
levels of gangliosides GMl, GDla, GTlb; c) GM3, GMl, 
and GM2 synthases expression. Therefore, the determin- 
ation of ganglioside pattern could significantly improve 
the prediction of clinical outcome and possibly help the 
design of the most appropriate therapeutic strategy. 

Conclusions 

In addition to be classified on the basis of characteristic 
genetic mutations, different types of melanomas could 
also be distinguished on the basis of ganglioside profile. 
This analysis that could be performed in histochemistry 
also on surgical specimens, could indicate the aggres- 
siveness of the tumor and, therefore, improve the predic- 
tion of the clinical outcome. 
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Additional file 1: Clinical data and ganglioside composition of 
melanoma cell lines and melanocytes. Origin of the cell line: VPG 
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